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1 Introduction 

The construction and study of classical (and quantum) dynamical r-matrices is currently 
undergoing extensive development. Various examples of such objects were recently discussed, 
for instance in ]l], 0, [|. However at this time one lacks a general classifying scheme such 
as exists in the case of constant classical r-matrices thanks to Belavin and Drin'feld H. A 
partial classification scheme has very recently been proposed || for dynamical r-matrices 
obeying the particular version of the dynamical Yang-Baxter equation || corresponding to 
Calogero-Moser models [0]. 

Consideration of such structures is thus particularly relevant to the study of Calogero- 
Moser || [| and relativistic Ruijsenaars-Schneider models jR], [UJ where they appear sys- 
tematically Their occurence and the particular form they assume, which we shall give in 
detail in this lecture, are due to the common nature of these models as hamiltonian reduc- 
tions of free or harmonic motions on particular symplectic manifolds: cotangent bundle of 
Lie algebras or Lie groups for rational or trigonometric Calogero-Moser models [T!|; double 



of Lie groups |13| or cotangent bundle to centrally extended loop group |14| for Ruijsenaars- 
Schneider models. The most general elliptic potentials are in turn associated to loop groups 



over elliptic curves [jLj, [13] and are crucial in understanding the algebraic resolution of these 
models 



15, 16 



This hamiltonian reduction procedure is also an important tool in the explicit construc- 
tion of classical r-matrices for 5C n -type systems, where the direct resolution of the intricate 
r-matrices equations, such as was done in [ I7| , becomes untractable. We shall exemplify 



such a construction in the case of trigonometric Calogero-Moser models |[Lq| . 

The plan of this lecture runs as follows. We first recall the essential results of clas- 
sical r-matrix theory and introduce the notations to be used throughout it. We then 
describe the construction of classical r-matrices for trigonometric Calogero-Moser models 
using the Hamiltonian reduction procedure. This gives us a general formula valid for all 
non-exceptional Lie groups. 

We finally give a systematic overall picture of the classical r-matrix structure of this type 
which were obtained by various authors for: 

a) elliptic, trigonometric and rational A n Calogero-Moser models [|HJ 



b) elliptic, trigonometric and rational A n spin Calogero-Moser models ]T9 



c) elliptic, trigonometric and rational A n Ruijsenaars-Schneider models. [20[ ^T], 



The case of spin Ruijsenaars-Schneider systems [11] was recently investigated in the 



rational case [[23|] but the general Hamiltonian structure still escapes understanding at this 
time. 
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To be complete we must indicate that two alternative approaches were recently described. 
One, using the same Lax matrices but leading to dynamical r-matrices of a different type 
(to be commented upon later) was developed in (see contribution by Pr. Braden). The 
other one uses conjugated Lax matrices which allow to eliminate the dynamical dependence 
in the r-matrix at the cost of introducing a more complicated Lax operator. It was developed 
directly in the quantum case using the formalism of intertwining vectors |2f| (see contribution 
by Pr. Hasegawa). 



2 Preliminaries 

First of all we need to recall four essential features of the classical r-matrix formalism (see 



2(| for a textbook presentation of the hamiltonian theory of classical integrable systems.). 

We consider a generic dynamical system described by a set of coordinates {xj} and mo- 
menta {pi}, i = 1 • • -n; a Poisson structure {} and a Hamiltonian h(xi,pj). 

1- Liouville theorem 

The existence of n algebraically independent, globally defined, Poisson commuting quan- 
tities such that the hamiltonian h belongs to the ring generated by this set, guarantees the 
existence of a canonical transformation (xi,pj) —>■ 9j) linearizing the equations of motion 
271] . Further assumptions on the topological structure of the phase space allow more precise 



statements on the geometrical interpretation of the transformed variables known as action- 
angle variables p^j. 

2- Lax pair formulation 

The Lax pair formulation of a dynamical system is the giving of two elements of a Lie 
algebra Q, L(x,p) and M(x,p) such that the equations of motion for x iy Pi are equivalent to 



the isospectral evolution (Lax equation) [29] 



f = [4M] (1) 

It follows that the adjoint-invariant quantities TrL n are time- independent. In order to 
implement Liouville theorem onto this set of possible action variables we need them to be 
Poisson-commuting. This is ensured by the classical r matrix structure. 

3- The r-matrix structure 

Defining the decomposition of the Lax operator on a basis {t a } of the Lie algebra Q as 
L = Y^a l a t a , the Poisson commutation of the ad-invariant TrL n is equivalent to the existence 
of an object r 12 (x,p) e Q\ <8> Q2 hereafter known as a classical r-matrix J3D], [JT], such 
that: 

{Li, L 2 } = £{Z°, l b }t a ®t b = [r 12 , Lx] - [r 21 , L 2 ] (2) 

a,b 

It must immediately be remarked that such an object is by no means unique. Moreover 
there is no one-to-one correspondance between a given dynamical system and the Lie algebra 
in which its Lax representation is defined; a same dynamical system may have several Lax 
representations and several r-matrix structures. 
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4- The classical Yang-Baxter equation 

The Poisson bracket structure (§) obeys a Jacobi identity which implies an algebraic con- 
straint for the r-matrix. Since r depends a priori on the dynamical variables this constraint 
takes a complicated form : 



[Li, [n 2 , r 23 ] + [ri 2 , r i3 ] + [r 32 , r i3 ] + {L 2 , r i3 } - {L 3 , r 12 }] + cycl. perm. = (3) 

Relevant particular cases of this very general identity are obtained when: 

a) r is independant of x, q. One is then lead to the general non-dynamical Yang Baxter 
equation |3(J [31], |32], p3| : 

[ri 2 , ^23] + [ru, ria] + [r 32 , r 13 ] = (4) 

b) If furthermore r is antisymmetric under permutation of the two copies of the algebra 
Q one obtains the better known and much studied [|3"U] , || form: 

[ri2, r 23 ] + [ri2, r i3 ] + [r 13 , r 23 ] = (5) 

c) If on the contrary r is dynamical, the supplementary terms in (|3|) may take a completely 
algebraic form due to the specific structure of the Lax operator. For instance the Calogero- 
Moser models lead to a self-contained algebraic equation for a classical r matrix depending 
only on position- type canonical variables x v : 



[ri2, r 13 ] + [n 2 , r 23 ] + [r 13 , r 23 ] - £ h^-^-r 23 + £ hf^-r 13 - ]T hf^-r 12 = (6) 

where the set {/i^} is a choice of basis for the Cartan algebra of £ acting on the representation 
space z. This equation was first derived in || and a classification scheme of its solutions was 
proposed in 0, closely connected to the general algebraic scheme in [f|. 

Note by the way that there exists canonical examples of dynamical r matrices, obtained 
by using the well-known higher Poisson bracket construction for any integrable dynamical 



system p4| , pq] starting from a constant r-matrix. In particular the quadratic Sklyanin 
bracket, where the Poisson structure of a Lax matrix becomes {Li, L 2 } = [7*12, L\ <S>£ 2 ], may 
also be described as a linear structure with a dynamical r matrix i?i 2 = ri 2 L 2 . We shall 
see on the example of the Ruijsenaars- Schneider system that dynamical "linear" r-matrix 
structures may also give rise to dynamical "quadratic" structures. The initial dynamical 
r-matrices themselves, however, are themselves not of the Sklyanin type. 



3 Hamiltonian reduction and r-matrices 
3.1 General hamiltonian reduction 

We begin by recalling some well-known facts concerning the Hamiltonian reduction of dy- 
namical systems whose phase space is a cotangent bundle [3B|. Let M be a manifold and 
N = T*M its cotangent bundle. N is equipped with the canonical 1-form a whose value at 
the point p G T*M is ir*p where 7r is the projection of N on M. If a Lie group G acts on 
M, each element X G Q (the Lie algebra of G) generates a vector field on M that we shall 
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denote X.m at the point m G M. It lifts to a vector field on N leaving a invariant. We shall 
also denote X.p G T p (N) the value at p G A" of this vector field, so that the Lie derivative 
Cx.p oi of the canonical 1-form vanishes. 

iV is a symplectic manifold equipped with the canonical 2-form to = —da To any function 
(or Hamiltonian) H on N we associate a vector field Xh such that dH = ix H uj and conversely 
since u is non-degenerate. 

The Hamiltonian associated to the above vector field X.p, X G Q reads: 

Hx (p) = «x. P a = a (Xp) (7) 

For any two functions F, G on N one defines the Poisson bracket {F, G} as a function 
on AT by: 

{F,G} = lu(X f ,X g ) (8) 

The Poisson bracket of the Hamiltonians associated to the group action has a simple expres- 
sion. In fact the group action is Poissonnian, i.e.: 

{Hx, Hy} = H[X,Y] (9) 

Obviously the application Ig?-* Hx{p), any foranyp G AT, is a linear map from to 
the scalars and so defines an element of Q* which is called the momentum at p G N. 

One then restricts oneself to the submanifold of N with fixed momentum /i i.e. such 
that: = V' 1 ^) 

Due to equation (0) and the invariance of a the action of the group G on N is transformed 
by V into the coadjoint action of G on (?* 

V(g.p)(X) = a(g.g- 1 Xg.p) = Ad* g V(p)(X) (10) 

where the coadjoint action on an element £ of (?* is defined as: 

Ad^(X)=£(^ 1 X( 7 ) 

The stabilizer G^ of /x G acts on A^. The reduced phase space is precisely obtained by 
taking the quotient (assumed well-behaved): 

= N./G, (11) 

It is known that this is a symplectic manifold. 

We then need to compute the Poisson bracket of functions on JF^. These functions are 
conveniently described as G^ invariant functions on A^. To compute their Poisson bracket 
we first extend them arbitrarily in the vicinity of N^. Two extensions differ by a function 
vanishing on N^. The difference of the Hamiltonian vector fields of two such extensions is 
controlled by the following: 

Lemma 3.1 Let f be a function defined in a vicinity of and vanishing on N^. Then the 
Hamiltonian vector field Xf associated to f is tangent to the orbit G.p at any point p G A^. 

As a consequence of this lemma we have a method to compute the reduced Poisson 
bracket. We take two functions defined on A^ and invariant under G M and extend them 
arbitrarily. Then we compute their Hamiltonian vector fields on N and project them on 
the tangent space to by adding a vector tangent to the orbit G.p. These projections are 
independent of the extensions and the reduced Poisson bracket is given by the value of the 
symplectic form on N acting on them. 
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Proposition 3.1 At each point p G one can choose a vector Vf.p G Q.p such that 
Xj + Vf.p G T p (N fl ) and Vf.p is determined up to a vector in Q^.p. 

One finally gets the consistent general formula for reduced Poisson brackets: 

Proposition 3.2 The reduced Poisson bracket of two functions on Tp, can be computed using 
any extensions f, g in the vicinity of according to: 

{/.^reduced = {/,</} + \ ((V g . P ).f - (Vf. P ).g) (12) 

This is equivalent to the Dirac bracket. 

3.2 The case N = T*G 

Let now M = G be a Lie group; one uses the left translations to identify N = T*G with 
Gxg*. 

to G T*(G) — > (<?,£) where u = L*-^ (13) 
The Poisson structure on N = T*G is easily seen to be: 

{£(X),£(Y)} = -a[X, Y}) ; (X),g} = -g X ; {g,g} = (14) 

Geodesies on the group G correspond to left translations of 1-parameter groups (the 
tangent vector is transported parallel to itself), therefore Stig^y) = 0. This is a Hamiltonian 
system whose Hamiltonian is: H = | (£,£)• where we have identified Q* and Q through the 
invariant Killing metric. 

Here H is bi-invariant, so one can reduce this dynamical system using Lie subgroups Hl 
and Hr of G of Lie algebras Hl and Hr, acting respectively on the left and on the right on 
T*G in order to obtain a non-trivial result. 

Using the coordinates (<?,£) on T*G this action reads: 

( (h L ,h R ), (0,0 ) ^ (h L gh- R \Ad* h J) 

We have written this action as a left action on T*G, in order to apply the formalism developed 
in Section 3.1. 

The moments are: 

V L (g,0 = P H >A.d;t; V R (g,0 = -P<H*Z; V=(V L ,V R ) (15) 

where we have introduced the projector on 7i* of forms in Q* induced by the restriction of 
these forms to TC. 

3.3 The Calogero-Moser models 

We follow here the derivation of [T^J. Let us consider an involutive automorphism a of a 
simple Lie group G and the subgroup H of its fixed points. Then H acts on the right on G 
defining a principal fiber bundle of total space G and base G/H, which is a global symmetric 
space. Moreover G acts on the left on G/H and in particular so does H itself. We shall 
consider the situation described in Section 3.2 when Hl = Hr = H. The Hamiltonian 
of the geodesic flow on G/H is invariant under the H action allowing to construct the 
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Hamiltonian reduction which under suitable choices of the momentum leads to the Calogero- 
Moser models. As a matter of fact since the phase space of the Calogero model is non-compact 
one has to start from a non-compact Lie group G and quotient it by a maximal compact 
subgroup H so that the symmetric space G/H is of the non-compact type. The derivative 
of a at the unit element of G is an involutive automorphism of Q also denoted a. Let us 
consider its eigenspaces TC and K, associated with the eigenvalues +1 and —1 respectively. 
Thus we have a decomposition: 

Q = H@K, (16) 

in which TC is the Lie algebra of H which acts by inner automorphisms on the vector space 
K, (hICh- 1 = K). 

Let A be a maximal commuting set of elements of K. It is called a Cartan algebra of the 
symmetric space G/H. It is known that every element in /C is conjugated to an element in 
A by an element of H. Moreover A can be extended to a maximal commutative subalgebra 
of Q by adding to it a suitably chosen abelian subalgebra B of TC. We shall use the radicial 
decomposition of Q under the abelian algebra A: 

G = AQ)B Re a (17) 

e a , oe* 

These decompositions of Q exponentiate to similar decompositions of G. First G = KH 
where K = exp(/C). Then A = exp(^4) is a maximal totally geodesic flat submanifold of 
G/H and any element of K can be written as k = hQh~ x with Q G A and h G H. It follows 
that any element of G can be written as g = h\Qhi with hi,h 2 G H. 

Of course this decomposition is non unique. This non-uniqueness is described in the 
following: 

Proposition 3.3 If g — h\Qh 2 = h , 1 Q , h , 2 we have: h[ = hidT 1 ^ , h' 2 = h dh 2 and Q' = 
hoQliQ 1 where d G exp(B) = B and h G H is a representative of an element of the Weyl 
group of the symmetric space. So if we fix Q = exp(g) such that q be in a fundamental Weyl 
chamber, the only ambiguity resides in the element d G B. 

The reduction to Calogero-Moser models is then obtained by an adequate choice of the 
momentum fi = (fi L , fi R ) such that V — fi. We take fi R = so that the isotropy group of 
the right component is Hr itself. 

The choice of the moment \i L is of course of crucial importance. It must be fixed so that: 

• its isotropy group is a maximal proper Lie subgroup of H, so that the phase space 
of the reduced system be of minimal dimension but non trivial. 

• In order to ensure the unicity of the decomposition introduced in the Proposition ( |3.3|) 
on Np we need: 

H^DB = {0} (18) 
We choose a complementary maximal isotropic subspace C so that 

n = n fl ®B®c (19) 

and x is a non-degenerate skew-symmetric bilinear form on B © C, hence dim£> = 
dimC. Notice that C is defined up to a symplectic transformation preserving B. 
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• The reduced phase space has dimension 2 dim A 

We now construct a section S of the bundle over so that one can write: 

N„ = H^SH (20) 

To construct this section we take a point Q in A and an L G £/* such that the point (Q, £) is 
in Np. In this subsection we shall for convenience identify Q and Q* under the Killing form 
assuming that G is semi-simple. Moreover since the automorphism a preserves the Killing 
form, Ti and K are orthogonal, and P-n* reduces to the orthogonal projection on Ti. Since 
fi = we have L G /C and one can write: 

L = p + L(e a -a(e a )) (21) 

e a , 

where p G A. From equation (0) one gets: 

^ L = Ph (p + E k (Q^Q- 1 - Qcxie^Q- 1 )^ 

Since Q = exp(g), q & A we have Qe Q Q _1 = exp(a(g))e Q , and similarly Qc^e^Q^ 1 = 
exp(-a(g))er(e a ) . 

Then the above equation becomes: 

H L = J2l a smha(q) (e a + cr(e a )) (22) 

a 

One can choose the momentum of the form: fi L = J2 a 9a(^a + o"(e a )) namely /i L has no 
component in B, where the g a are such that is of maximal dimension (we shall see that 
it essentially fixes them, and obviously if g a ^ for any a equation (|IED is automatically 
satisfied) and we have shown the: 

Proposition 3.4 The couples (Q,L) with Q = exp(g) and 

„ smha(g) 

with p,q G A form a submanifold in of dimension 2 dirndl. 

Notice that L is just the Lax operator of the Calogero model and that the section S 
depends of 2 dim A parameters in an immersive way Hence one can identify iV M with the set 
of orbits of S under if M x H i.e. the set of points (g = hxQh 2 , £ = h 2 l Lh 2 ) with h\ G and 
h 2 G H uniquely defined due to condition fllTf). The variables p and q appearing in Q and L 
are the dynamical variables of the Calogero model and form a pair of canonically conjugate 
variables. 

We then compute the Poisson bracket of the functions on whose expressions on the 
section S are L(X) and L(Y) for X, Y G K. These functions have uniquely defined x H 
invariant extensions to given respectively by: 

Fx(g, =< f , h 2 1 Xh 2 >, F Y (g, =< f , /ij 1 Yh 2 > where 5 = h x Qh 2 

Notice that h 2 is a well-defined function of g in iV M due to condition (|1|). According to the 
prescription given in the section 3.1 we choose extensions of these functions in the vicinity 
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of N^. We define these extensions at the point p = (<?,£) G T*G by the same formulae in 
which h 2 is chosen to be a function depending only on g and reducing to the above-defined 
h 2 when p G iV M . Because of the non-uniqueness of the decomposition g = hiQh 2 outside 
Nfj, one cannot assert that the functions Fx, Fy are invariant under the action of H x H 
and we must appeal to the general procedure to compute the reduced Poisson brackets. 



The complete derivation with all its technical subtleties can be found in [18]. The final 



result gives the general r-matrix for trigonometric Calogero-Moser models in so-called dual 



form pOfl : 



Theorem 3.1 There exists a linear mapping R : K, — > 7i such that: 

{ L(X) , L(Y) } reduced = L ( [X, RY] + [RX, Y] ) (23) 

and R is given by: 

R(X) = V„fro (XI 4- 

where: 



R{X) = V 9 h 2 {X)+ l -D Q {V x ) (24) 



Proposition 3.5 On the section S with Q = exp(g) G A we have: 
For X G K i.e. X = X + Y,X a (e a - <re a ), X e A 

V g h 2 (X) = -h (X) + J2X a coth(a(q))(e a + ae a ) (25) 

a 

Here h Q (X) is a linear function from Q to B which is fixed by the condition: 

X L = h (X) - sin *2(q) (6a + a6a) eH ^® C ( 26 ) 
3.4 Two examples 

To illustrate the power of this method we now give two examples of r-matrices. The A n case 
had already been treated in W7\ and serves as a check on the validity of the derivation. The 
case of SU (n, n) Calogero-Moser model proved to be too intricate for a direct computation; 
however this method immediately gives its r-matrix. 

The standard Calogero-Moser model (SL(n)) is obtained by starting from the non 
compact group G = SL(n, C) and its maximal compact subgroup H = SU(ri) as first shown 
by fli~2"f . We choose the momentum \i L as described in Section (3.1) so that the isotropy 
group Hp be a maximal proper Lie subgroup of H. Obviously one can take fii of the form: 

Hl = * {w + - 1) (27) 

where v is a vector in C n such that v + v = 1, hence \xl is a traceless antihermitian matrix. 
Then g^^g^ 1 = fi^ if and only if gv = cv where c is a complex number of modulus 1. Hence 
= S(U(n — 1) x £7(1)) which has the above-stated property. 

In this case the automorphism a is given by a (g) = (g + )^ 1 (notice that we consider only 
the real Lie group structure), B is the group of diagonal matrices of determinant 1 with 
pure phases on the diagonal and A is the group of real diagonal matrices with determinant 



I. The property (18) is then satisfied as soon as the vector v has no zero component. As a 
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matter of fact, v is further constrained by \xl being a value of the moment map. Considering 
equation (|22"D we see that hl has no diagonal element, which implies that all the components 
of v are pure phases Vj = exp(i9j). These extra phases which will appear in the Lax matrix 
can however be conjugated out by the adjoint action of a constant matrix diag (exp(idj)) 
hence we shall from now on set Vj = 1 for all j. This is the solution first considered by 
Olshanetskii and Perelomov. 

The Lax matrix L is then given by Proposition (|3.4p and therefore 



P + ■ w 1 \( iE ki - iEi k ) 

^ smh(g fc - qi) 



The r-matrix can now be deduced straightforwardly from Proposition ( |3.5| ), after recon- 
verting the dual form where R is an endomorphism of the Lie algebra into the more usual 
direct form where R lives in the tensor product of the Lie algebra by itself. One ends up 
with: 



11 1 

R12 = coth (<?fc - qi)E u ®E lk + ~Y^ — -(E kk 1) <g> (E M - E lk ) 

This gives back the already known r-matrix of the Calogero model for the potential 1/ sinh(x), 
and the other potentials 1/ sin(x) and 1/x have similar r-matrices obtained by analytic con- 
tinuation. 

The SU (n, n) Calogero model is obtained by starting from the non-compact group 
G = SU(n,n). This is the subgroup of SL(2n,C) which leaves invariant the sesquilinear 
quadratic form defined by 

Q((u 1 ,v 1 ),(u 2 ,v 2 )) = (ut vt)j( U v 2 J=u+v 2 + vtu 2 (28) 

where Ui, Vi are vectors in C" and J is the matrix 



J 



1 

1 



The Lie algebra of SU (n, n) therefore consists of block matrices 

Q = {( a c J) \a = -d + , Tr (a + d) = 0, b + = -b, c + = -c] (29) 

where a, b, c, d are n x n complex matrices. 

We consider again the automorphism o : cr(g) = (g + )~ 1 , which can be consistently 
restricted to SU(n,n). Its fixed points at the Lie algebra level consist of block matrices 

H = {(° c C \ \a + = -a, Tr(a) = 0, c+ = -c} (30) 

This Lie algebra is isomorphic to the Lie algebra of S(U(n) x U(n)), the two w(n)'s being 
realized respectively by a + c and a — c. 
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The subalgebra B consists of matrices of the form (^) with c = 0, and a is a diagonal 
matrix of zero trace and purely imaginary coefficients. The Abelian subalgebra A consists 
of matrices of the form fl2"9|) with b = c = and a = —d being a real diagonal matrix. 

To perform the reduction, we choose as above p, R = and 

H L = i(vv + - 1) + i7 J (31) 

The vector v has again 2n components all equal to 1. 

Notice that in equation ([IT]) the parameter 7 is an arbitrary real number. This will lead 
to existence of a second coupling constant in the corresponding Calogero model. 

Then, from proposition (|3.4|), the Lax matrix is found to be: 



L = P + Yl ■ , / 1 r ( 1 - o) (iE i:j + iE j+ni+n ) 

+ g sinh(g! + 9j ) (1 " a)iiEi ' j+n + iEj ' i+n) 

+(7 + i)E-^(i-)(^) 

where p is a generic element of A of the form diagp^, — diagp^. 
The r-matrix is then computed straightforwardly: 

= 9 X! cot ' 1 (q , fc — qi)(E k i + E fe+rij ; +ri ) ® (i?^ — Ei +riik+n ) 

+ \ C0tll (gfc + qi)(E k> i +n + E k+n ,l) <8> (^+n,fc - -E«,fc+n) 

11 1 

+ "X! ■ , / AE kk + E k+riyk+n 1) ® (i^H — E k+n j +n ) 

£ kjkl smn^fc — qi) n 

+ - Y, . -, ; r(-^fcfc + E k+n k+n 1) <g) (E k i +n — E k+n i) 

2 ^ smh{q k + qi) n 

These dynamical r-matrices depend only on the dynamical variable g. The different 



approach advocated in p3 leads to r-matrices depending on both p and q variables, but on 



a smaller set of algebra generators. 



4 The dynamical r matrices of Calogero and Ruijse- 
naars models 

Dynamical r-matrices have been derived for the Calogero- Moser and (relativistic) Ruijsenaars 
Schneider models using either the technique described here or a direct method starting from 
an ansatz of the same form. We will now describe the results achieved in this way for 
A n models, and indicate interesting and sometimes deep connections between these various 
r-matrices. 

Let us start with Calogero-Moser models. The rational and trigonometric matrices 
were described in the previous section.. The elliptic case was solved by Sklyanin and by 



10 



Braden et al. [O. The Lax matrix reads: 



TV N 



i=l 



Here one has set: 



Z(x,A) 



cr(x + A) 



V{x) = p(x) 



a(x) er(A) ' 

where cr and p are Weierstrass elliptic functions. The classical r-matrix reads: 



(32) 



(33) 



A' 



N 



ris(K AO = J2 ^ - A 4 ) e *j ® + 2 XI %i> A 4 ) ( e « + e ii) ® 



V 



-[C(A - //) + C(A*)] 2 e i{ ® e«. 



8=1 



(34) 



Note that a spectral parameter is now present in L and r. This particular formulation 



of the spinless elliptic case is due to Krichever ||15||. The other known Lax formulation due 



to Olshanetskii and Perelomov [12] has no spectral parameter but requires a p and q de- 



pendance in the r-matrix, which was only recently given in |24j and has a totally different 
algebraic form. 

The spin Calogero-Moser models were introduced in ||. The Lax operator for the 
elliptic case reads: 



N 



N 



L W =J2Pi e H + A) fij e { 



(35) 



8=1 



where fij are spinlike variables with the Kirillov-Poisson bracket structure: 



{fij, fkl} — 2 i^il fjk + Ski fij + fijk fil + $lj fki)- 



(36) 



One then needs to introduce a parametrization of so as to be on a coadjoint orbit of 
SU(N): Introducing vectors 

(&)i=l...JV With & = (C")a=l-r 

(r)i)i=i-N with 7ft = (^ a ) a= i... r 



with the Poisson brackets 



= 0, K, = 0, {&, I,}} = 5, 



06 1 



we parametrize as follows 



a=l 



The phase space now becomes a true symplectic manifold. 



(37) 



(3f 
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The Hamiltonian takes the form 

N j N 

F= nE^-nE k U v (Qa), Qij = ( li- Qj ( 39 ) 

Z i=l Z i,j=l 

with the Weierstrass function as elliptic potential, as in the spinless case. 
The classical r-matrix then reads: 



I N I N 

^ »,J=1 ^ i,J=l 

1 ^ 

- r [C(A + fi) + C(A - //)] £ e - ® e «- (40) 
^ i=i 

Trigonometric and rational cases can be derived from the elliptic case by taking suitable 
limits pi. The spinless case can also be derived from the spin case by taking r = 1 and 



introducing a further Hamiltonian reduction by the action of U(l) as a phase on the vectors 
£i,T]i. The supplementary terms in the spinless r matrix arise from the conjugation of the 
Lax matrix required to bring it in canonical shape ( |32"D after elimination of the vector-like 
degrees of freedom. Let us finally remark that this r-matrix structure has yielded a number 
of important developments: exact classical Yangian symmetry |H| (a quantum version of 



it had been found beforehand, using heavy direct algebraic computations [STfl); quantum 
version of the dynamical r-matrix using the shifted version of the quantum Yang-Baxter 
equation described in |39], [3£| [4(J . 



The spinless relativistic RS models are described by the hamiltonian: 

N 

H = mc 2 Y, (coshfl,-) J] f(q k - q 3 ) (41) 



where 



i 9 2 Y /2 
f(q) = I 1 + — I (rational) 

/ a 2 \ 1/2 
fW = ( 1 + 2 (hyperbolic) 



x sinh , 

f(q) = (X + uV(q)) (elliptic), V = Weierstrass function (42) 

Here the canonical variables are a set of rapidities i — 1 • • • iV} and conjugate positions 
qi such that {0i, qj} = dij. 

The dynamical system admits a Lax representation with the Lax operator: 



exp (36 j ■ C jk (qj - q k ) 




1/2 



n/(a 



(43) 
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where {e^} is the usual basis for N x N matrices; / was given in (52) and 

7 

Cjkiq) = (rational) 

7 + iq 

Cjk{q) = ( cosh— q + ia sinh— q J (trigonometric) (44) 



2" 2 
® (q + J,\ ) 
$( 7 , A) 



C M = 1Z. [\ ( elli P tic ) ( 45 ) 



Again the function $ is defined as: 

+ £ (46) 
a{x)a{X) 

where a is the Weierstrass function. 

The elliptic r-matrix structure is better written as a quadratic expression in terms of the 
Lax operator [22|: 



{Lx(X), L 2 (fj)} = (la (A) ® L 2 (//))ai(A, //) - a 2 (A, /i)(Li(A) <g> L 2 (aO) 
+(l<g>L 2 (/x))s 1 (A,//)(L 1 (A)®l) 

-(L 1 (A)®l)s 2 (A,/i)(l®L 2 (/i)) (47) 

Here one defines: 



at(X, ft) — a(X, fi) + w , Si(A,//) = s(A, fi) — w 
a 2 (X, fi) = a(X, ft) + s(X) — s*(fi) — w , s 2 (A, fi) = s* (ft) + w (48) 

The matrices a and s are obtained from the r-matrix of the elliptic Calogero-Moser model 
given in fl34]) as r(X, fi) = a(X, fi) + s(X) where a is the skew-symmetric matrix: 

N 

a(X, fi) = -C(A -fi)J2 E kk ®E kk ~Y, ®(qj ~ Qk, A - ® ^ (49) 
fc=i fc^j 

and s, s* are non-skew-symmetric matrices independent of the second spectral parameter: 

N 

s(X, ft) = c(A) XI ^ ® E kk + X ~ Ik, X)E jk ® £ fcfe 
fe=i fc^j 

s*(A,/i) = C(A)X)- B fcfc®" E7 ** + S $ (3j" - Q , *' A )" E; *fc® £ y* ( 50 ) 
fc=i fc^i 

and finally w is a supplementary matrix, independent of the spectral parameters: 

w = SC(?fc- Qj)Ekk® Ejj (51) 
Mi 

This r-matrix structure is a Sklyanin-type bracket (although realized in the more generic 
case of an initially dynamical r-matrix) obtained from the Calogero-Moser r-matrix structure 

13 



viewed as a linear bracket. This can be interpreted from the fact that RS models are 
obtained not only as hamiltonian reductions from current algebras on elliptic curves but also 
alternatively |13[ as hamiltonian reductions from Heisenberg double of Lie groups . In 
this case the initial Poisson structure on the large phase space is itself a quadratic bracket 
instead of the canonical initial linear (Kirillov) bracket which is the natural structure on the 
cotangent bundle of a Lie group. This relation is maintained throughout the hamiltonian 
reduction procedure and the final r-matrix structures are essentially connected in the same 
way. 

The previously obtained r-matrices [^TJ can be obtained from this one by sending one 
period of the elliptic functions to infinity and suitably conjugating the Lax pair in such a way 
as to get a completely symmetric expression in terms of the momenta 0j. On the other hand 
the r-matrix found in |20] cannot be easily inserted in this scheme. In fact it correspond to 
a very specific value of the parameters where the Lax matrix becomes completely symmetric 
and tne r-matrix may only then take this very special form. 

The classical r-matrices admit a quantization scheme on the same lines as the Calogero- 
Moser case J7J. It was developed in [4"2"|. 



Finally a word about the spin RS dynamical system. They were introduced in [11]. It 



is not clear at this time how to define a consistent hamiltonian structure in the most general 
case although the rational case was solved recently [p^ |. The task is indeed easier here since 
there exists a duality symmetry j|3| connecting the rational RS model to the trigonometric 
CM model for which the spin model is well known. Let us finally mention that there exists 
a general scheme to obtain Hamiltonian structures from Lax representations using the tools 
of algebraic geometry |0J (see also Pr. Krichever's contribution to this colloquium) and 
this scheme now appears to be the most promising way to get these elusive Hamiltonian 
structures. 
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